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反 向 蒙特 卡 罗 方法 在 卫星 辐射 
分 析 中 的 研究 与 应 用 
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摘 要 反 向 蒙特 卡 罗 方法 (AMC/RMC) 是 Geant 4 中 一 个 强 有 力 的 偏 置 技术 . 粒子 产生 于 包含 灵敏 体 的 反 向 
界面 并 在 几何 体 中 被 反 向 追踪 直至 外 界面 或 溢出 能 量 闵 值 ， 其 计算 时 间 只 集中 于 对 结果 有 贡献 的 粒子 径 迹 .与 正 
向 蒙特 卡 罗 方法 (FMC) 相 比 , RMC 更 高 效 , 特别 是 当 灵 人 敏 体 远 小 于 几何 体 其 他 部 分 和 外 界面 时 ,其 优势 更 明显 . 
通过 RMC 与 EMC 的 比较 , 验证 了 RMC 应 用 于 卫星 辐射 剂量 分 析 的 准确 性 . RMC 与 SHIELDOSE2 和 SSAT 
的 比较 说 明了 RMC 是 高 精度 卫星 辐射 剂量 的 优选 方法 . 
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Abstract The Reverse Monte Carlo (RMC) method is a powerful biasing technique available in 
Geant 4, also known as the Adjoint Monte Carlo method. In this method, particles are generated 
on the external boundary of the sensitive part of the geometry and then tracked backward in the 
geometry till they reach the external source surface, or exceed an energy threshold. By this way 
the computing time is focused only on particle tracks that are contributing to the tallies. The RMC 
method is much rapider than the Forward Monte Carlo (FMC) method when the sensitive part of the 
geometry is smaller compared to the rest of the geometry as well as the external source. Compared 
with the FMC method, accuracy of the RMC method applied to the radiation dose of the satellite 
is verified. Comparisons of RMC with the SHIELDOSE2, SSAT show that RMC is the preferred 
method of high-precision satellite radiation dose. 
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0 引言 


空间 辐射 环境 -3 的 主要 来 源 包括 地 球 辐射 
带 、 太 阳 高 能 粒子 、 银 河 宇宙 射线 等 . 由 于 这 些 辐射 
源 的 存在 , 空间 辐射 环境 十 分 复杂 . 辐射 带 粒 子 和 太 
阳 高 能 粒子 等 不 仅 局 限 在 材料 表面 发 生 相互 作用 , 能 
量 高 的 粒子 会 穿 过 材料 表面 进入 材料 内 部 , 影响 航天 
器 内 部 载荷 的 正常 运行 . 研究 分 析 辐 射 效应 对 航天 带 
的 辐射 防护 具有 重要 意义 . 

前 向 ( 正 向 ) 蒙特 卡 罗 方 法 (EMC)I 起 源 于 20 
世纪 40 年 代 , 是 一 种 随机 抽样 技巧 , 粒子 直接 从 外 部 
射 入 并 正 向 追踪 . 为 得 到 期 望 的 精确 模拟 结果 , 蒙特 
卡 罗 模 拟 通 常 需要 大 量 的 计算 时 间 .， 当 灵敏 体 与 剩 
余 几 何 体 以 及 外 界面 相 比 均 很 小 时 , 很 多 计算 时 间 花 
费 在 对 探测 器 信号 并 无 贡献 的 粒子 群 上 . 这 种 情况 下 
反 向 蒙特 卡 罗 方 法 (RMC/AMC) 能 用 于 加 速 计算 . 
反 向 蒙特 卡 罗 方 法 最 早 被 成 功 应 用 是 在 20 世纪 60 
年 代 , 核 科学 家 通过 反 向 蒙特 卡 罗 方 法 得 到 核反应 堆 
某 些 部 分 辐射 剂量 的 快速 评估 [入 引 ， 在 空间 系统 辐 
射 分 析 方面 反 向 蒙特 卡 罗 方 法 还 被 用 于 NOVICE 代 
码 中 9. 

常见 的 蒙特 卡 罗 方 法 仿真 软件 主要 有 MORSE， 
EGS*，MNCP*，NOVICE***，Geant 4** 等 ， 其 
中 NOVICE 和 Geant4 可 实现 反 向 蒙特 卡 罗 模 拟 . 
NOVICE 商用 软件 目前 无 法 获取 , Geant 4 给 出 了 反 
向 蒙特 卡 罗 方 法 的 开源 , 但 只 提供 两 种 简单 能 谱 的 反 
向 蒙特 卡 罗 辐 射 计算 中 .在 中 国 反 向 蒙特 卡 罗 方 法 
主要 用 于 医学 .飞行 器 、 机 械 制 造 等 红外 热 辐 射 特性 
研究 , Dong, Qi, Shuai 等 & 2 在 这 方面 有 较 多 研究 
成 果 , 但 在 航天 器 剂量 分 析 领 域 的 应 用 相对 较 少 . 本 
文 基于 Geant 4 的 RMC 模块 实现 了 任意 能 谱 的 反 
向 蒙特 卡 罗 剂 量 分析 , 通过 与 FMC, SHIELDOSE2， 
SSAT 等 进行 比较 , 研究 反 向 蒙特 卡 罗 方法 用 于 高 精 
度 剂 量 计 算 的 可 行 性 . 


1 基本 原理 与 方法 


1.1 _ Geant 4 RMC 模块 
反 向 蒙特 卡 罗 方 法 (Reverse/Adjoint Monte 


*http://rcwww.kek.jp/research/egs/ 
**https://mcnp.lanl.gov/ 
***http://www.empc.com/novice.php 
****http://geant4.cern.ch/ 
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Carlo, RMC/AMC) 是 Geant 4 中 一 个 强 有 力 的 偏 置 
技术 . 粒子 在 几何 体 灵 人 敏 部 分 的 反 向 界面 ( 见 图 1 界 
面 4) 上 产生 , 之 后 在 几何 体 中 被 反 向 追踪 直至 到 达 
外 界面 或 者 溢出 能 量 闵 值 ， 这 样 计算 时 间 只 和 针对 于 
对 结果 有 贡献 的 粒子 径 迹 .， 当 几何 体 的 灵敏 部 分 相 
比 其 他 部 分 和 外 界面 ( 见 图 1 界面 B) (扩散 的 且 不 像 
激光 束 ) 很 小 时 , RMC 方法 比 正 向 蒙特 卡 罗 方法 快 
很 多 . 目前 RMC 方法 在 Geant4 中 只 用 于 实现 部 分 
电磁 过 程 , 即 电子 、 质 子 和 离子 的 反 向 离散 电离 ( 通 
过 电子 离 化 和 韧 致 辐射 以 及 质子 和 离子 电离 得 到 的 
持续 能 量 获得 ) 以 及 电子 的 反 向 离散 韧 致 辐射 、 反 向 
光电 效应 、 反 向 康 善 顿 散射 和 近似 多 重 散射 ， RMC 
是 Geant 4 中 一 个 重要 模块 , 从 版 本 9.3 加 入 , 这 里 
所 使 用 的 版 本 为 Geant 4.9.5. 

1.2 反 向 蒙特 卡 罗 方 法 

如 图 1 所 示 , 反 向 粒子 ( 反 向 电子 、 反 向 7 粒 
子 等 ) 在 反 向 界面 4 持续 产生 , 带 有 随机 的 位 置 、 能 
量 (1/B 分 布 ) 和 方向 ， 反 向 界面 为 用 户 定义 体 或 
球 的 外 界面 , 包含 一 个 或 几 个 灵敏 体 , 相对 于 整个 几 
何 体 来 说 是 小 的 .用户 可 设 定 反 向 界面 的 最 大 和 最 小 
值 ， 反 向 初级 粒子 产生 后 在 几何 体 中 被 反 向 追踪 直 
至 到 达 一 个 用 户 定义 的 外 界面 B (球面 或 立体 界面 ) 
或 在 达到 用 户 定义 的 最 高 能 量 限制 前 被 终止 , 这 个 能 
量 限 制 代表 外 界面 的 最 大 能 量 . 在 反 向 踪迹 中 , 反 向 
过 程 发 生 在 反 向 粒子 在 几何 体内 被 追踪 或 散射 或 转 
换 成 其 他 类 型 的 反 向 粒子 处 ， 一 个 反 向 追踪 阶段 对 
应 一 个 Geant 4 事件 . 

如 图 1 中 路 径 1 所 示 , 当 一 个 反 向 粒子 达到 外 
界面 4 时 , 其 权重 、 类 型 、 位 置 和 方向 被 记录 , 同时 
一 个 与 最 后 产生 的 初级 反 向 粒子 类 型 相同 、 能 量 相 
同 但 方向 相反 的 常规 粒子 产生 , 在 灵敏 区 域 被 正 向 追 
踪 , 像 正 向 蒙特 卡 罗 方 法 模拟 一 样 进行 权重 抽样 . 一 
个 正 向 追踪 阶段 对 应 一 个 Geant 4 事件 . 

产生 于 反 向 界面 并 在 模拟 中 考虑 的 反 向 和 正 向 
粒子 类 型 表 是 一 个 在 物理 表 中 声明 过 的 反 向 过 程 函 
数 . 例如 如 果 只 考虑 电子 和 7 电磁 过 程 , 则 只 将 反 
向 电子 和 反 向 7 粒子 作为 初级 粒子 . 这 样 一 个 反 向 
事件 将 被 分 为 4 个 G4 事件 , 即 反 向 电子 的 反 向 追 
踪 、 正 向 电子 的 正 向 追踪 、 反 向 7 粒子 的 反 向 追踪 
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和 正 向 7 粒子 的 正 向 追踪 . 因此 100 个 反 向 事件 需 
考虑 400 个 G4 事件 .如 果 将 质子 的 离 化 也 考虑 进 
来 , 反 向 质子 和 正 向 质子 也 作为 初级 粒子 , 则 100 个 
反 向 事件 需 考虑 600 个 G4 事件 . 

1.3 ”真实 环境 模拟 处 理 

Geant4 中 只 能 实现 两 种 能 谱 分 布 , 即 简单 的 指 
数 分 布 和 须 指 数 分 布 ， 而 空间 环境 却 是 复杂 多 变 的 ， 
为 实现 广泛 应 用 这 里 将 针对 外 部 环境 真实 能 谱 进行 
分 析 . 

外 界 环境 采用 的 是 辐射 带 模 型 AE8D3 在 太阳 
峰值 时 处 于 低 轨 道 (近地点 400km, 远地点 500 km) 
高 度 的 捕获 电子 微分 能 谱 , 如 图 2 所 示 . 在 反 向 蒙特 
卡 罗 模 拟 中 需 对 图 2 所 示 的 任意 能 谱 采 用 线性 插值 
方法 , 插值 多 项 式 


p(E) = JEo) 十 


f(B1) — fl(Eo) 
Bi — Eo 


(B—Eo). (1) 


图 1 反 向 蒙特 卡 罗 模块 在 Geant 4 中 的 实现 过 程 
Fig.1 Schematic view of an adjoint/reverse 


simulation in Geant 4 
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图 2 外 部 环境 卫星 轨道 能 谱 (AE8Max) 
Fig.2 Satellite particle energy spectrum (AE8Max) 
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1.4” 归 一 化 处 理 

正 向 蒙特 卡 罗 方 法 得 到 的 模拟 结果 是 针对 每 个 
事件 数 给 出 的 , 其 对 应 于 归 一 化 到 外 部 源 发 出 的 一 个 
粒子 积分 通 量 .实际 剂量 需 进 一 步 通过 公式 处 理 得 
到 真实 粒子 数 才能 计算 分 析 . 

归 一 化 粒子 数 NN 主要 与 能 谱 分 布 和 角 分 布 有 
关 . 假设 给 定 外 部 积分 能 谱 下 (单位 cm 2.s sr 下)， 
则 Nm = 64m?R*， 这 里 给 定 的 是 离散 全 向 微分 能 
谱 65' (单位 cm ?.s-1), 有 


= [Ben Ei)(fi 


firi(Biti — Ei)|, 


fi+1)1 


即 6’ = 了 于 (已 + — Ei)(fit fiti). 
由 此 得 出 正 向 归 一 化 粒子 数 


Ni = $nR. (2) 


与 正 向 蒙特 卡 罗 方 法 模拟 类 似 , 反 向 蒙特 卡 罗 
方法 模拟 通过 抽样 微分 截面 计算 反 向 权重 因子 归 一 
化 到 真实 情况 .在 反 向 蒙特 卡 罗 方法 中 ,对 每 一 个 
散射 初级 粒子 和 有 反 向 反应 的 产物 均 考虑 一 个 相应 的 
独立 反应 . 一 个 给 定 的 反 向 反应 的 抽样 反 向 微分 
截面 O20o1()/08009 等 于 带 有 倒置 的 发 射 体 和 产 
物 (或 散射 发 射 体 ) 能 量 的 正 向 微分 截面 , 即 

O20t(E) 

OF'0N 
其 中 , O20(B')/9E609 为 正 向 反应 的 微分 截面 , 2 给 
出 了 产物 或 散射 发 射 体 相 对 于 发 射 体 的 方向 ， 不 考 
虑 其 他 偏 置 和 权重 修正 , 反 向 权重 为 


O20(E’) 


E00 局 


品 z 0) 
DPI101 
oe-Je® Di Bor)dryy,. (4) 


W =Waelo DhlBis)dr Nya dxdEodfo. 


式 中 ,Wa 与 探测 表面 的 反 向 初级 粒子 的 抽样 有 关 ， 
Wh 为 外 部 通 量 的 归 一 化 因子 . 外 部 通 量 的 归 一 化 是 
通过 乘 以 反 向 粒子 到 达 外 界面 的 权重 来 实现 的 , 同时 
伴随 相应 的 方向 通 量 . 因此 到 达 外 界面 的 反 向 权重 


W =dE1df1dS4 cos Gae-)o! Dii(Bur)de Ny. 
B20 (Bo) 
B10n 


dzxdEodfoe do" Sie(Bor)dr f( po, 10). 
(5) 


206 
2 辐射 剂量 分 析 


2.1 RMC 与 FMC 电子 屏蔽 分 析 
2.1.1 ”几何 模型 与 环境 

所 模拟 的 结构 体 为 半径 10cm 屏蔽 厚度 1mm 
的 简单 球体 ， 灵 敏 体 为 位 于 球 心 的 半径 5mm 高 
度 5mm 的 圆柱 体 , 垂直 方向 距离 球 心 6cm 和 一 6cm 
处 各 有 一 块 长 宽 均 为 4cm 厚度 0.25mm 的 Ta 板 . 精 
度 作为 截断 误差 取 为 5%. 计算 机 CPU 型 号 为 Intel 
Corei5-2400, 主 频 3.10 GHz. 

精度 定义 为 正 向 和 反 向 蒙特 卡 罗 模 拟 时 计算 沉 
积 在 灵敏 体 的 能 量 所 设 定 的 预期 精确 度 . 即 当 计算 出 
的 相对 统计 误差 沉积 能 量 达到 该 精度 时 运行 中 止 并 
记录 结果 ， 否 则 继续 运行 按照 用 户 指定 的 事件 数 处 
理 . 在 默认 情况 下 预期 精度 被 设置 为 0, 这 意味 着 运 
行将 不 会 中 止 . 
2.1.2 ”模拟 结果 对 比分 析 

将 反 向 蒙特 卡 罗 横 拟 结果 与 正 向 蒙特 卡 罗 作 对 
比 , 得 到 的 总 能 量 分 别 为 132.72MeV 和 133.86 MeV， 
即 剂 量 为 1.16041x10-6rad 和 1.17035x10-6rad， 
将 正 向 蒙特 卡 罗 结 果 作 为 参考 值 , 相对 误差 为 0.35%. 
由 于 Geant4 本 身 反 向 物理 过 程 处 理 带 来 的 误差 及 
其 他 误差 , 正 向 与 反 向 蒙特 卡 罗 模 拟 对 比 误 差 达到 
约 10% 都 是 可 信服 的 . 因此 反 向 蒙特 卡 罗 方 法 与 正 
向 蒙特 卡 罗 方 法 在 电子 屏蔽 剂量 分 析 中 的 准确 性 是 
可 比拟 的 . 

图 3~5 分 别 为 给 定 几 何 体 在 捕获 电子 能 谱 作用 
下 正 向 和 反 向 蒙特 卡 罗 得 到 的 模拟 结果 对 比 . 图 3 为 
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图 3 灵敏 体 沉积 能 量 与 初始 电子 能 量 的 变化 关系 


Fig.3 Energy deposited in the sensitive volume versus 


primary electron energy 
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灵敏 体 中 沉积 能 量 随 初始 电子 能 量 的 变化 , 正 向 和 反 
向 蒙特 卡 罗 结 果 符 合 较 好 . 图 4 为 进入 灵敏 体 的 电子 
通 鞭 随 能 量 的 变化 , 图 5 为 进入 灵敏 体 的 7 通 量 随 
能 量 的 变化 . 

误差 大 的 部 分 原因 在 于 初始 能 谱 高 能 电子 较 少 
导入 外 部 真实 环境 的 模拟 离散 插值 也 会 带 来 一 部 分 
误差 .在 极 少数 情况 下 ,到达 外 部 源 时 反 向 轨迹 可 
能 会 产生 一 个 很 高 的 权重 . 虽然 这 不 显著 发 生 , 但 往 
往 会 影响 模拟 结果 ， 反 向 多 重 散 射 和 反 向 币 至 辐射 
在 Geant4 中 的 处 理 也 有 一 定 误差 . 对 于 反 向 多 次 散 
射 , 使 用 与 正 向 相同 的 模型 . 这 种 近似 使 得 反 向 与 正 
向 蒙特 卡 罗 方 法 模拟 之 间 存 在 差异 . 对 于 件 致 辐射， 
反 向 与 正 向 蒙特 卡 罗 模 拟 模 型 中 使 用 的 微分 截面 之 
间 的 区 别 是 在 反 向 蒙特 卡 罗 模 拟 中 > 100 keV 的 产 
生 界 面 有 更 高 的 7 通 量 . 反 向 过 程 应 利用 方向 微分 
截面 对 反 向 次 级 粒子 进行 抽样 并 计算 反 向 微分 截面 
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图 4 灵敏 体 中 电子 通 量 随 能 量 的 变化 


Fig.4 Fluence of electron entering the sensitive volume 
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图 5 灵敏 体 中 7 通 量 随 能 量 的 变化 


Fig.5 Fluence of 7 entering the sensitive volume 
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而 事实 上 使 用 的 微分 截面 是 通过 数值 推导 由 正 向 模 
型 提供 的 方向 截面 在 次 级 7 能 量 上 获得 的 , 正 反 疝 
微分 截面 不 匹配 . 

图 6 给 出 的 是 电子 离散 能 谱 环 境 精 度 取 5% 屏 
蔽 厚度 d 取 0.1mm 情况 下 灵敏 体 大 小 变化 影响 正 
向 与 反 向 蒙特 卡 罗 方 法 的 模拟 结果 圆柱 体 灵敏 体 
大 小 由 灵敏 体 底面 半径 和 高 度 决定 . 

图 6 图 (a) 为 两 种 方法 针对 不 同 灵敏 体 大 小 对 
沉积 剂量 的 影响 . 由 图 6(a) 可 知 , 灵敏 体 大 小 对 沉积 
剂量 的 影响 很 大 , 尺寸 越 小 剂量 沉积 越 多 , 相对 而 言 
计算 难度 越 大 ， 而 正 向 和 反 向 两 种 蒙特 卡 罗 方法 模 
拟 结果 符合 较 好 , 误差 较 小 , 均 在 6% 以 内 , 这 说 明 自 
定义 轨道 能 谱 应 用 于 电子 屏蔽 的 Geant 4 工具 包 蒙 
特 卡 罗 计 算是 可 行 的 , 通过 与 正 向 蒙特 卡 罗 模 拟 结果 
的 对 比 说 明 反 向 蒙特 卡 罗 用 于 电子 屏蔽 剂量 计算 准 
确 度 是 可 信服 的 . 

图 6(b) 为 两 种 方法 对 比 模拟 的 相对 误差 与 灵敏 
体 大 小 之 间 的 关系 . 由 于 Geant4 本 身 反 向 物理 过 程 
处 理 带 来 的 误差 及 其 他 误差 使 得 正 向 和 反 向 蒙特 卡 
罗 模 拟 对 比 误 差 在 约 10% 内 都 是 可 信服 的 . 图 6(b) 
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中 可 见 误差 浮动 区 间 为 0~6%, 平均 误差 3.69%, 结 
果 较 好 . 

图 6(c) 给 出 了 两 种 方法 在 精度 取 5% 而 屏蔽 厚 
度 d 取 0.1mm 情况 下 要 达到 预期 精度 的 计算 时 间 随 
灵敏 体 大 小 的 变化 . 对 于 正 向 蒙特 卡 罗 方 法 而 言 , 灵 
敏 体 尺 寸 越 小 耗 时 越 入 ， 从 几 分 钟 变化 到 几 十 小 时 ， 
跨越 3 个 数量 级 , 而 反 向 蒙特 卡 罗 方 法 基本 耗 时 不 
变 , 稳定 在 几 秒 钟 . 随 着 屏蔽 厚度 和 精度 的 变化 , 反 向 
蒙特 卡 罗 方 法 计算 时 间 也 较 稳定 而 正 问 蒙特 卡 罗 方 
法 则 不 稳定 . 这 是 反 向 蒙特 卡 罗 方 向 的 一 个 优点 . 在 
模拟 中 还 发 现 精度 与 计算 时 间 成 负 相关 , 精度 越 小 ， 
计算 时 间 越 长 , 正 向 蒙特 卡 罗 方法 耗损 的 时 间 越 多 . 
这 可 用 于 电子 的 屏蔽 计算 ， 质 子 模拟 过 程 计算 时 间 
较 短 , 这 是 因为 质子 追踪 路 径 多 为 直线 ,而 电子 散射 
较 多 . 

用 于 鉴定 运算 时 间 的 品质 因子 灵敏 值 FOM ( 定 
义 符号 SFom) 为 

i 7 (6) 

其 中 , 为 计算 沉积 能 量 的 精确 度 , 了 为 以 分 钟 为 单 
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图 6 灵敏 体 大 小 对 沉积 剂量 和 计算 时 间 的 影响 


Fig.6 Different sizes of the sensitive volume affect on computed dose and computing time 
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位 的 计算 所 消耗 的 时 间 ， 本 文 精度 R 取 5%, 正 向 
蒙特 卡 罗 方 法 和 反 向 蒙特 卡 罗 方 法 的 计算 时 间 分 别 
为 41085.4s 和 19.609s， 品 质 因子 分 别 为 0.58 和 
1224， 这 意味 着 要 达到 相同 的 精度 反 向 蒙特 卡 罗 方 
法 比 正 向 蒙特 卡 罗 方 法 快 大 约 2110 倍 . 这 便 是 反 向 
蒙特 卡 罗 方 法 的 最 大 优势 之 处 . 

由 此 可 知 , 反 向 蒙特 卡 罗 方 法 是 计算 剂量 精度 较 
高 且 高 效 的 方法 . 
2.2 RMC 与 SHIELDOSE2 和 扇 段 分 析 法 比较 

目前 用 于 计算 剂量 的 常用 方法 是 简单 的 球 模 
型 和 平板 模型 法 、 较 复 杂 的 扇 段 法 (射线 追踪 
法 ) 以 及 三 维 蒙特 卡 罗 方 法 . 简单 的 球 模型 和 平板 
模型 法 是 将 卫星 结构 简化 为 简单 的 几何 结构 ， 计 
算 简单 几何 屏蔽 结构 中 的 辐射 剂量 ， 代 表 工 具 例 
如 SHIELDOSE203- 芍 ， 户 段 分 析 法 (射线 追踪 法 ) 
是 考虑 卫星 三 维 屏蔽 结构 , 但 简化 粒子 与 材料 的 相互 
作用 , 代表 工具 例如 SSATH5. 三 维 蒙特 卡 罗 方 法 既 
考虑 粒子 与 材料 的 相互 作用 过 程 , 又 可 以 处 理 复杂 的 
三 维 屏蔽 结构 , 所 以 精度 更 高 , 而 且 反 向 蒙特 卡 罗 比 
正 向 蒙特 卡 罗 更 高 效 . 
2.2.1 与 SHIELDOSE2 比较 

SHIELDOSE2 是 航天 工程 中 常用 的 辐射 剂量 计 
算 软 件 ， 可 以 计算 四 种 简单 屏蔽 下 的 辐射 剂量 ， 包 
括 有 限 平 板 、 半 无 限 平板 、 实 心 球 和 球 壳 屏 蔽 . 选 
择 实心 球 模型 与 RMC 进行 比较 . 选择 实心 球 材 料 
为 Al (密度 2.7g'cm 3), 目标 灵敏 球 的 材料 为 Si ( 密 
度 2.33 gcm-3). 由 于 Geant 4 中 不 人 允许 载荷 有 重 丢 ， 
在 反 向 蒙特 卡 罗 方 法 模拟 中 需要 挖 空 比 灵敏 球 大 一 
点 的 区 域 . 选择 挖 空 半 径 为 灵敏 体 球 半 径 . 外 界 环境 
采用 的 是 电子 辐射 带 模型 AE8 在 太阳 峰值 时 处 于 地 
球 同 步 轨道 高 度 的 粒子 能 谱 数据 , 时 间 周 期 为 一 年 . 

图 7 给 出 的 是 SHELDOSE2 和 反 向 蒙特 卡 
罗 方 法 模拟 得 到 的 捕获 电子 沉积 剂量 随 Al 屏蔽 实 
心 球 半 径 变 化 的 曲线 .Al 屏蔽 球 尺寸 由 屏蔽 半径 
决定 . 在 同步 轨道 上 质子 虽然 对 总 剂量 也 有 贡献 ， 
但 是 贡献 最 大 的 是 电子 , 而 且 其 随 半径 变化 最 显著 ， 
所 以 对 比 的 主要 是 捕获 电子 沉积 剂量 . 可 以 看 出 ， 
二 者 随 Al 屏蔽 球 尺寸 变化 趋势 基本 相同 , 半径 越 
小 则 屏蔽 厚度 越 小 沉积 剂量 越 大 ， 随 半径 增 大 变 
化 趋势 越剧 烈 . 例如 当 Al 屏蔽 球 半径 取 3mm 灵 
敏 体 球 半径 取 lnm 时 ，SHIELDOSE2 捕获 电子 沉 
积 剂量 为 3.27800x104rad,，RMC 计算 的 出 的 剂量 
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为 32462.69rad, 相对 误差 为 0.98%. 当 Al 屏蔽 球 半 
径 取 7mm 灵敏 体 球 半径 取 1nm 时 , SHIELDOSE2 
捕获 电子 沉积 剂量 为 3.639 00x102rad, RMC 计算 的 
剂量 为 219.67rad, 相对 误差 为 65.65%. 将 RMC 计 
算 结果 作为 参考 值 , 相对 误差 变化 分 析 如 图 8 所 示 . 

图 8 表明 了 SHIELDOSE2 和 RMC 两 种 方法 模 
拟 得 出 的 沉积 剂量 之 前 的 相对 误差 与 Al 屏蔽 球 半径 
变化 之 间 的 关系 . 当 Al 屏蔽 球 半径 范围 在 0.1~5mm 
内 时 两 者 误差 较 小 , 随 着 屏蔽 球 半 径 增 大 , 误差 在 半 
径 为 3mm 时 达到 一 个 最 小 峰值 , 半径 5mm 之 后 继 
续 增 大 则 误差 也 随 之 增 大 ， 在 半径 较 小 或 较 大 时 误 
差 都 较 大 , 误差 从 0.98% 到 65% 变化 不 等 , 平均 误差 
为 30.52%. 

前 面 通过 与 正 向 蒙特 卡 罗 方 法 的 比较 已 经 说 明 
反 向 蒙特 卡 罗 方法 在 电子 屏蔽 计算 剂量 时 的 准确 度 
是 可 信 的 . SHIELDOSE2 是 利用 预先 计算 的 电子 和 
质子 蒙特 卡 罗 模 拟 结果 推算 出 来 的 剂量 na9, 图 8 表 
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明 在 屏蔽 尺寸 较 小 或 较 大 时 ，SHIELDOSE2 的 计算 
误差 会 增加 . 因此 , 反 向 蒙特 卡 罗 方 法 是 高 精度 计算 
剂量 时 的 更 优选 择 . 
2.2.2 与 扇 段 分 析 法 比较 

扇 段 分 析 法 "是 航天 器 设计 中 计算 辐射 剂量 
的 常用 方法 ,这 种 方法 简化 了 粒子 与 材料 的 相互 作 
用 , 将 所 关心 位 置 的 全 向 立体 角 划 分 为 若干 小 的 记 
面 , 在 每 个 书面 内 发 出 者 干 条 射线 , 直线 追踪 射线 路 
径 的 物质 , 并 积分 路 径 上 的 等 效 厚度 , 对 等 效 厚 度 求 
平均 得 到 该 扇面 的 屏蔽 厚度 , 利用 该 扇面 的 等 效 屏蔽 
厚度 在 剂量 深度 曲线 中 播 值得 到 该 扇面 的 剂量 贡献 . 

SSATU5] 为 欧 空 局 (ESA) 基于 Geant4 开发 
的 户 段 屏蔽 分 析 工 具 ， 目 前 已 得 到 广泛 应 用 . 本 文 
将 RMC 的 结果 与 SSAT 的 结果 进行 了 比较 . 

RMC 三 维 结构 体 模型 如 图 9 所 示 ， 选 择 长 
度 20 cm 屏蔽 厚度 2mm 的 正方 体 壳 作为 卫星 载体 ， 
半径 1mm 的 球形 载荷 作为 灵敏 体 , 材料 均 为 Al. 卫 
星 有 效 载荷 4 个 , 垂直 方向 距离 中 心 士 6 cm 处 各 有 
一 块 长 宽 均 为 5cm 厚度 1mm 的 平板 , 材料 为 Ta, 水 
平 沿 工 轴 方 向 距离 中 心 土 5 cm 处 各 有 一 个 长 度 2 cm 
的 正方 体 盒子 , 材料 为 Al. 截断 误差 精度 取 5%. 

SSAT 几何 模型 与 RMC 基本 相同 , 唯一 不 同 是 


图 9 RMC 结构 体 模型 
Fig.9 Satellite Structure of RMC 
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少 了 灵敏 体 ，SSAT 的 屏蔽 厚度 计算 方法 有 法 向 传 
输 (NORM) 和 直线 传输 (SLANT) 两 种 , 这 里 采用 的 
是 SLANT 方法 . RMC 外 界 环境 采用 的 是 电子 辐射 
带 模 型 AE8 在 太阳 峰值 时 处 于 地 球 同步 轨道 高 度 时 
间 周 期 为 一 年 的 离散 能 谱 . SSAT 采用 的 是 如 图 7 所 
示 SHIELDOSE2 的 实心 球 模型 剂量 深度 曲线 . 

不 同 灵敏 体位 置 对 沉积 剂量 的 影响 列 于 表 1, 平 
均 误 差 为 31.28%. 位 置 a 是 卫星 结构 中 心 (0, 0, 0)， 
以 RMC 方法 计算 结果 作为 参考 值 RMC 与 SSAT 
相对 误差 为 49.10%, 误差 较 大 .位置 5 是 (0, 一 35， 
0) mm， 上 距离 了 轴 负 半 轴 盒子 中 心 15mm, RMC 
与 SSAT 相对 误差 为 28.46%. 位 置 c 为 (0，0， 
50) mm, 距离 2 轴 正 半 轴 Ta 板 中 心 10mm, 两 者 
相对 误差 为 24.82%. 位 置 4 为 (0, 0, 90)mm, 位 于 2 
轴 卫 星 Al 立方 体 党 和 Ta 板 之 间 , 距离 Al 立方 体 
党 10mm, 距离 Ta 板 30 mm, RMC 与 SSAT 相对 误 
差 为 22.74%. 


3 结论 与 讨论 


通过 自 定义 能 谱 环 境 的 电子 屏蔽 和 质子 屏蔽 分 
析 , 比较 RMC 与 FMC 两 种 方法 的 沉积 剂量 和 计算 
时 间 . 

(1) 自 定 义 能 谱 计 算得 到 的 灵敏 体 中 沉积 能 量 随 
初始 电子 能 量 的 变化 , RMC 和 FMC 结果 符合 较 好 . 
以 FMC 结果 作为 参考 值 电子 屏蔽 计算 相对 误差 
为 0.85%. 

(2) RMC 和 FMC 两 种 蒙特 卡 罗 方 法 模拟 结果 
曲线 符合 较 好 , 误差 较 小 . 不 同 灵敏 体 大 小 两 种 方法 
的 相对 误差 浮动 区 间 为 0~6%, 平均 误差 3.69%, 结 
果 较 好 . 

(3) 对 比分 析 两 种 方法 的 计算 时 间 , 电子 屏蔽 计 
算 时 RMC 比 FMC 速度 提高 2110 倍 , 质子 屏蔽 计 


表 1 不 同位 置 RMC 和 SSAT 剂量 对 比 
Table 1 Comparison of radiation dose between RMC and SSAT method at different locations 


位 置 点 SSAT 相对 误 
沉积 能 量 /MeV 剂量 /rad 剂量 /rad 差 /(%) 

a 2.496 54 x10-8 353 18.82 526 60 49.10 

b 2.253 08 x10-8 318 74.56 409 47.6 28.46 

C 2.31871 x10-8 328 03.04 409 45.7 24.82 

d 2.433 17 x10-8 34422.31 422 51.3 22.74 
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算 时 速度 提高 43 倍 . 在 相同 情况 下 , RMC 和 了 FMC 
的 计算 时 间 随 灵敏 体 大 小 的 变化 而 变化 . FMC 耗 时 
久 且 变化 大 , 而 RMC 耗 时 基本 不 变 . 随 着 屏蔽 厚度 
和 精度 的 变化 反 向 蒙特 卡 罗 方 法 计算 时 间 较 稳定 而 
正 向 蒙特 卡 罗 方 法 则 不 稳定 

与 FMC 模拟 结果 的 对 比 表 明 RMC 自 定义 
能 谱 运 用 于 卫星 辐射 剂量 分 析 是 可 行 且 可 信 的 . 而 
且 RMC 在 计算 时 间 上 有 很 大 优势 , 不 仅 高 效 而 且 
稳定 . 

对 比 RMC 与 SHIELDOSE2 和 SSAT 的 沉积 剂 
量 , 得 到 以 下 结果 . 

(1) RMC 和 SHIELDOSE2 两 种 方法 的 剂量 深度 
曲线 基本 相符 . 随 Al 屏蔽 球 尺 寸 变化 趋势 基本 相同 
半径 越 小 则 屏蔽 厚度 越 小 , 沉积 剂量 越 大 . 以 RMC 
计算 结果 作为 参考 值 ， 相 对 误差 随 半径 变化 有 浮动 ， 
中 间 小 两 头 大 , 平均 误差 为 30.52%. 

(2) RMC 和 SSAT 对 比 使 用 自 定义 卫星 模型 , 考 
虑 处 于 几何 体 不 同位 置 的 灵敏 体 , 并 计算 分 析 沉 积 剂 
量 . 以 RMC 结果 为 参考 值 , SSAT 与 RMC 相对 误 
差 平均 为 31.28%. 

SHIELDOSE2 实现 简单 屏蔽 球 的 剂量 计算 ， 
SSAT 忽略 了 粒子 与 材料 的 相互 作用 . 与 这 两 种 方法 
计算 结果 的 比较 可 知 , 反 向 蒙特 卡 罗 方 法 是 可 信和 且 精 
度 较 高 的 计算 方法 . 

反 向 蒙特 卡 罗 方 法 是 一 种 快速 且 高 精度 的 计算 
剂量 方法 . 由 于 Geant 4 发 展 的 局 限 性 , 后 期 将 研究 
重合 卫 星 结构 和 改进 反 向 蒙特 卡 罗 模 拟 的 粒子 和 有 反 
向 物理 过 程 . 
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